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Abstract 

An outline is given how to formulate a relativistic unitarized constituent quark model 
of mesons in momentum space, employing harmonic quark confinement. As a first step, 
the momentum-space harmonic-oscillator potential is solved in a relativistically covariant, 
three-dimensional quasipotential framework for scalar particles, using the spline technique. 
Then, an illustrative toy model with the same dynamical equations but now one qq and one 
meson-meson channel, coupled to one another through quark exchange describing the 3 Pq 
mechanism, is solved in closed form on a spline basis. Conclusions are presented on how to 
generalize the latter to a realistic multichannel quark/meson model. 



1 INTRODUCTION 



Mesonic qq states are the simplest color-singlet configurations in QCD. However, physical mesons 
appear to be much more complicated objects, in view of the quite disparate spectra of especially the 
light scalars and pseudoscalars, and the large variety of states ranging from OZI-stable quarkonia 
to extremely broad, often badly established resonances. Because of these complexities, theoretical 
approaches usually focus on parts of the meson data only. Thus, one has quarkonium potential 
models (see e.g. Ref. Jl|), relativistic constituent quark models (see e.g. Ref. @), chiral quark 
models (see e.g. Ref. 0), heavy-quark effective theory [||], the linear a model (LaM) ||, Nambu- 
Jona-Lasinio models ||, (unitarized) chiral perturbation theory |7[]), meson-exchange models (see 
e.g. Ref. M), and so on. While it lies outside the scope of this short paper to describe in detail 
the respective merits and drawbacks of all these approaches, on can generally state that a good 
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reproduction of either spectra or scattering data is achieved, often at the cost of a non-negligible 
number of free fit parameters, but without accomplishing a unified description of mesonic spectra 
and meson-meson scattering. In our view, however, spectra and scattering of mesons are inexorably 
intertwined concepts, considering that most mesons are resonances, decaying strongly into final 
states of two or more lighter mesons. 

Therefore, we believe one should treat the mechanisms responsible for quark confinement and for 
mesonic decay on an equal footing. This we have realized already in an essentially non-relativistic 
(NR) framework, employing the coupled-channel Schrodinger formalism, with several confined qq 
channels and many free two-meson channels. For a detailed discussion of the configuration-space 
model, we refer to a separate contribution |§ to this workshop, so let us just emphasize here some 
important non-perturbative features. From the spectroscopic point of view, meson-loop effects 
on the qq spectra are included to all (ladder) orders, even those resulting from closed meson- 
meson channels. Conversely, the strong influence of s-channel qq states on non-exotic meson- 
meson scattering is taken into account. Furthermore, bound states, resonances, and scattering 
observables all follow directly from an analytic S-matrix derived in closed form, thus dispensing 
with any kind of perturbative methods, which could be highly untrustworthy for the often very 
large couplings involved. As a matter of fact, perhaps the most remarkable conquest of the model is 
its parameter-free prediction of the light scalar mesons ||10|| , including the k(800) (K$) which 
has very recently been experimentally confirmed P2]| . Precisely for the light scalars, perturbation 
theory manifestly breaks down [[|, p~3|j . 

Other coupled-channel, also called unitarized, approaches are due to Eichten et al. (Ref. |L|, 
first paper), Bicudo & Ribeiro (Ref. ||, fourth paper), and Tornqvist with collaborators (Ref. jl4l 
and references therein). The conclusions of the latter two models are in qualitative agreement 
with ours, in the sense that large shifts are found owing to unitarization, even for (bare) states 
below the lowest decay threshold. However, in the case of the scalar mesons, important differences 
between our model and especially the one of Ref. |T3| come to light, as the k(800) is not found in 
the latter reference (see Ref. for a detailed comparison). 

Notwithstanding its successes, our coordinate-space model has, of course, a number of short- 
comings. First of all, the light pseudoscalar mesons are quite badly off, except for the kaon. This 
may be due to the largely NR formalism without Dirac spin structure , the disregard of chiral 
symmetry, and possibly also the neglect of scalar-pseudoscalar two-meson channels. Moreover, for 
most processes, the experimental meson-meson phase shifts are only roughly and not very accu- 
rately reproduced, which is no wonder in view of the neglect of final-state interactions from direct 
t-channel meson exchange, and the absence of any "fitting freedom". Then there are problems 
with pseudothresholds in heavy-light systems, owing to the non-covariance of minimally relativis- 
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tic equations in configuration space. Finally, too many, viz. 3, parameters — albeit a very modest 
number as compared to other approaches — are needed for the 3 Po transition potential, due to 
the lack of a microscopic description in terms of quark exchange. 

Nevertheless, we are convinced one can deal with these shortcomings by reformulating the 
model in momentum space, in a covariant three-dimensional (3D) quasipotential framework. Here, 
one could object why we do not attempt to tackle the full Bethe-Salpeter (BS) equation. The 
reason is twofold: first of all, the numerical effort to solve the resulting coupled 4D integro- 
differential equations with many channels would be enormous; but perhaps more importantly, the 
confinement mechanism employing potentials is an inherently 3D concept, possibly leading to a 
pathological, non-confining behavior when generalized to four dimensions (see e.g. Ref. f2|, third 
paper). Anyhow, a covariant formulation would also allow to include, at a later stage, the Dirac 
spin structure of the quarks, and address the issue of dynamical chiral-symmetry breaking and 
generation of constituent mass. The resulting multichannel equations in momentum space we shall 
show to be explicitly solvable, provided we stick to harmonic confinement and work on a spline 
basis. 

This paper is organized as follows. In Sec. 2 the good old harmonic oscillator (HO) is solved in 
momentum space using the spline technique, for two different 3D relativistic equations. In Sec. 3 a 
two-channel toy model, with one harmonically confined qq channel coupled via quark exchange to 
one free two-meson channel, is formally solved in closed form, employing again a spline expansion. 
Conclusions and an outlook on how to further develop the model are presented in Sec. 4. 



In this section, we review the well-known 3D HO, first in coordinate space, and then in momentum 
space. In the latter representation, a 3D relativistic generalization is straightforward. 

2.1 Non- Relativistic Harmonic Oscillator 

The equal-mass radial Schrodinger equation (SE) for a local potential reads 



2 RELATIVISTIC HARMONIC OSCILLATOR 




(i) 



If we take the HO potential 



define x = 




(2) 




(3) 
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The spectrum of the 3D HO potential is very well known, viz. (/ = 0) 

E = tu(2n+^), n = 0,l,2,..., (4) 

leading to the eigenvalues 3, 7, 11, . . . for Eq. (3). 

Alternatively, we can work in momentum space. The homogeneous Lippmann-Schwinger (LS) 
equation for the bound-state vertex function reads 

m = J V(P,P') Goif; E) T{p') , (5) 
where the free two-body Green's function is given by 

G \p') = E- P — + ie. (6) 

m 



Now define the wave function ^ = G T. Then Eq. (5) becomes 

d 3 p' 

W) 



Go\m(P) = J ^V(p,p')*(?) , (7) 
or equivalently 

This is just the SE in momentum space. In this representation, the HO potential is given by 



V(p,p') = ~ptu; 2 S^(p-p')A f: (9) 



which is local as in coordinate space. Substituting into (8) yields 

r 2 

m 2 



/, _ - ^aA*(p) = E*(p). (10) 



Taking again I = and introducing v(p) = p\l/(p), with p = \p\, we obtain 



p 2 1 2 rf 2 

T 

m 2 dp 



v (p) = Ev &) • ( n ) 



If we define the dimensionless variable y = p/ ' JJEj and write m = 2/i, we get 



{y*-^}"M = ^- (12 > 

This is identical to Eq. (3). Note that now the roles of the kinetic and the potential term are 
interchanged. 
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2.2 Blankenbecler— Sugar— Logunov— Tavkhelidze (BSLT) Equation 



The momentum-space equal-mass (scalar) BSLT [|15|] equation for the wave function has exactly 
the same form as in the LS case of Eq. (7), but now with center-of-mass (CM) Green's function 



(13) 



where E v = \Jp 2 + m 2 , and s is the total CM energy squared s = P = (p\ + p?) = (2m + E) . 
Note that G^ SLT can be derived in a completely covariant fashion, and written down explicitly for 
an arbitrary frame [16|]. In the NR limit we have E p « m and s ~ Am 2 + AmE, recovering the LS 
Green's function (6). With the new Gq 1 , we get instead of Eq. (8) 

:rl "■ ■ d 3 p' _, s s — 4m 2 



p " m 



v(p,p')W') 



*(p). 



In terms of the above-defined dimensionless variable the BSLT equation becomes 



(14) 



1 



d 2 



y 



to 

2m' 



,2 dy 2 



v{y) 



Am 2 



2muo 



v{y) ■ 



(15) 



From this equation we readily see that deviations from the LS case are determined by the parameter 
uj/2m. Furthermore, in the limit u/2m — > 0, the eigenvalues 3, 7, 11, . . . get recovered. However, 
there are additional relativistic effects due to s ^ 4m 2 + AmE. 

2.3 Equal-Time (ET) Equation 



The ET, also called Mandelzweig- Wallace |17j, Green's function can be obtained from the BSLT 
(or Salpeter) equation by adding an extra term Gq to the two-body Green's function Gq which 
approximates crossed-ladder contributions. This extra term becomes exact in the eikonal limit, or 
when one of the particles gets infinitely heavy. Therefore, the covariant ET equation manifestly 
has the correct one-body limit, i.e., one recovers either the Klein-Gordon or the Dirac equation, 
for scalar or spin-|ll/2 particles, respectively. In the scalar case, the CM ET Green's function is 
given by 



m c 



2E 2 - - 

v 4 



(16) 



The final form of the S'-wave ET equation for the HO potential in momentum space reads 

s — Am 2 



y 2 - 



l 



2 - 



1 



1 

— UT 

2m 



y 



2mu 



1 - 



2m 



1 + — 

2m 



d /) 



(17) 
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2.4 Non-Zero Orbital Angular Momentum 

For I 7^ 0, the dimensionless LS equation becomes 



d 2 1(1 + 1) 2 \ , s E i ^ 

+ V ' + x 2 \ u(x) = — u{x) . (U 



dx 2 < :2 



x* | -u; 



In order to numerically solve the LS and corresponding BSLT, ET equations, we substitute u(x) 
x l u(x). Then 

d 2 21 ^ d x 2 1 _, , 



) , , + ^(l-x^) + x 2 \u(x) = —u(x) . (19) 
[ ax z x z dx J ^ 

The function u(x) has the right boundary conditions to allow an easy numerical solution using an 

expansion on a spline basis (for details, see Ref. [18] and second paper of Ref. ||). Thus, the NR 



HO spectrum is recovered with great accuracy, viz. 

E = uo (2n + £ + -), n,/ = 0,1,2,... . (20) 

In the case of the BSLT and ET equations, the same can be done by replacing everywhere 

-— - -M- -) i- m\ 

dy 2 ~^ y 2 { y dy' dy 2 ' { ' 

The LS, BSLT, and ET results for I = are depicted in Fig. [I]. We see that the relativistic 
corrections are relatively modest for the BSLT equation, but can become huge in the ET case, 
especially for excited states. The latter large effect has to do with the sign change of the ET 
propagator (16) as a function of the relative momentum, and has also been observed in the Dirac 
case (see Ref. 0, second paper). However, this overshooting of the ET equation will strongly 
depend on the chosen Dirac structure of the confining force in a realistic model, and hence should 
not cause any worries at this stage. 

3 UNITARIZED MODEL IN MOMENTUM SPACE 

Consider now a simple two-channel model of one confined qq system coupled to one free meson- 
meson channel. In a unitarized picture, this can formally be described by a 2 x 2 T matrix, with 
the dynamical equation (LS, BSLT, or ET) 

/ T 00 T m \ ( Koo Vbi \ / Voo Voi \( G W T 00 T 01 \ 
{ T w T n J ~ ^ V 10 V u J + { V w V u ) { G 1 ){ T w T n ) ' 

Equation (22) is diagrammatically represented in Fig. |2]. Dashed lines are mesons, solid lines with 
arrows are (anti)quarks, and curly exchanges symbolize confinement. Furthermore, open circles 
represent meson-quark-antiquark vertices, and dots stand for three-meson vertices. Note that the 
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2E/uj 




6- 



4- 



' 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 

uj/m q 

Figure 1: Dimensionless energy levels of the HO for the LS, BSLT, and ET equations, as a function 
of uj/m q ; first three radial levels are depicted (I = 0). 
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Tn = Vn + VnGiTn + VioGqToi 



Figure 2: Graphical representation of Eq. (22); see text for explanation. 



communication between the qq and two-meson channels takes place through quark exchange. This 
will give rise to a generally non-local 3 P transition potential in configuration space, as soon as the 
quarks and mesons have different masses. The precise form of such a potential will also depend on 
possible form factors to be attributed to the various vertices. 

In order to solve Eq. (22), we must realize that no asymptotic scattering of quarks is possible, 
so Too, T i, and T w must in the end be eliminated from the equations: only Tn can be on shell. 
If we now take again harmonic confinement, Eq. (22) becomes a set of coupled integro-differential 
equations. Omitting for clarity angular-momentum indices, we get e.g. for T i 

1 r k 2 dk 

T i(p,p') = V 01 (p,p') - -^ 2 A p Go(p)T 01 (p,p') + J — V 01 (p,k)G 1 (k)T 11 (k,p > ) . (23) 

Note that p' is an essentially dummy variable in the equations, so we can choose a fixed value for 
it, say p. If we now expand T 01 and T n on a spline basis and discretize the initial momentum p, 
then T i can be expressed in terms of T u as follows. 

N 

T i(Vi,v) = E44(pJ , < = V i(Pi,p); (24) 

N 

T n (k,p) = E4^ii(k); (25) 

4\ = {l + i^ 2 A Pi G (p 4 )}4i(P,) ; (26) 

At = J^Vo^p^G^sUk). (27) 
Combining Eqs. (24)-(27), we can write 

A i • t i = v i + A n • t n t 01 = A^ 1 ■ {v i + A n • t n } , (28) 

where the boldface capitals and small letters are obvious shorthands for matrices and vectors, 
respectively. Turning next to T n , we get the integral equation 

T n (p,p) = Vn(p,p) + /^Vn(p,k)Gi(k)T 11 (k,p) + |^\/ 1 o(p,k)Go(k)T 01 (k,p). (29) 
Similarly to Eqs. (24)-(27), we define now 

B l A = s{M) - j^V 11 (p i ,k)G 1 (k)s{ 1 (k) ; (30) 
r k 2 r/k 

= J ^y 10 ( Pi ,k)Go(k) Soi (k), v\, = MPi,P)- (31) 
Then we can express tn in terms of itself, using Eq. (28): 
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Bn • tn = vn + B i • t i = vn + B i • Aqi 1 • {v 01 + A n • t n } . (32) 

(28) 

This finally leads to 

{B n - Boi • Aqi 1 • An} • tn = v u + B i ■ A^ 1 ■ v i , (33) 
t n = {B n - Boi • A^ 1 • An}" 1 ■ {v n + B i • A^ 1 • v i} . (34) 

Thus, we have obtained a closed-form expression for the T matrix (hence the S matrix), con- 
taining the complete information of the confinement-plus-scattering process. Needless to say that 
all these algebraic manipulations only make sense if the involved integrals actually exist and the 
spline expansions converge. Yet, we are confident this will be the case when appropriate boundary 
conditions are chosen, since the same spline expansion we are using here has been successfully 
applied to few-body scattering calculations in momentum space [ i"9fl . Nevertheless, in the Dirac 



case, to be dealt with in the future, form factors will have to be included for the meson-quark- 
antiquark vertices so as to make the integrals involving three quark propagators convergent. In 
any case, such form factors are physically meaningful, and should therefore already be considered 
in the approximation with spinless fermions, whenever doing phenomenology. 

Bound states and resonances are given by zeroes in the determinant of the inverted-matrix factor 
right after the equal sign in Eq. (34). Resonance poles in the second Riemann sheet can be searched 
for by analytic continuation to complex energies, using e.g. contour-rotation methods. Moreover, 
from the structure of the equations one readily sees that the inclusion of final-state interactions, 
due to t-channel meson exchange in the meson-meson channel, simply amounts to taking a non- 
vanishing Vn in Eq. (29), which will not lead to any additional numerical effort. Finally, the 
generalization to many scattering channels and also more qq channels is straightforward, at the 
expense of bigger matrices. 



4 CONCLUSIONS AND OUTLOOK 

In the foregoing, we have demonstrated that the formulation of a unitarized quark/meson model 
with harmonic confinement is relatively easy in momentum space. Thus, the 3 Po mechanism 
is dynamically described through P-wave quark exchange. Furthermore, the inclusion of final- 
state interactions via t-channel meson exchange is straightforward, and does not complicate the 
structure of the equations. The numerical resolution of these equations using splines looks very 
convenient, which should result in (approximate) analytic solutions, thereby facilitating an analytic 
continuation to complex energies in order to search for resonance poles. 
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Our first priority is, of course, to achieve a total control of the numerics, for real as well as com- 
plex energy. Then a simple toy model will be worked out so as to reproduce known coordinate-space 
results. Furthermore, cases problematic in the configuation-space formulation will be studied, such 
as heavy-light mesons. Also, the influence of t-channel meson exchange on meson-meson phase 
shifts will be investigated. In a somewhat more distant perspective, we intend to do a lot of phe- 
nomenology, and use the feedback to further refine the model. 
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